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BOUNDARY -LAYER  MEASUREMENTS  IN  3.84-  BY  10-INCH 
SUPERSONIC  CHANNEL 
By  Paul  F.  Brinich 


SUMMARY 

Boundary-layer  measurements  -were  made  in  the  transonic  and 
supersonic  regions  of  a  channel  having  maximum  cross-sectional 
dimensions  3.84  hy  10  inches  and  designed  hy  potential-flow  meth¬ 
ods  for  a  uniform  Mach  number  of  2.08  in  the  test  section.  At 
inlet  pressures  from  37  to  13  inches  of  mercury  absolute ,  turbu¬ 
lent  boundary  layers  were  observed  throughout  the  channel;  at  an 
inlet  pressure  of  5  inches,  laminar  boundary  layers  were  observed 
near  the  channel  entrance  with  turbulent  layers  downstream. 

A  comparison  of  the  experimental  and  theoretically  computed 
boundary  layers  at  the  high  inlet  pressures  showed  good  agreement 
when  empirical  friction  coefficients  were  evaluated  from  Reynolds 
numbers  based  on  the  kinematic  viscosity  of  the  air  at  the  wall. 
Despite  this  agreement  between  experiment  and  theory,  local  dif¬ 
ferences  in  rates  of  boundary -layer  growth  still  existed  that  are 
attributed  to  secondary  flows  in  the  boundary  layer. 

At  low  inlet  pressures,  substantial  increases  in  the  boundary- 
layer  rates  of  growth  with  an  uneven  development  of  the  boundary 
layer  along  the  bottom  wall  of  the  channel  were  observed. 

.  Experimental  and  empirical  skin-friction  coefficients  were  in 
poor  agreement  at  all  inlet  pressures .  Secondary  flows  in  the 
boundary  layer  caused  by  static-pressure  gradients  transverse  to 
the  stream  direction  are  believed  to  be  the  reasons  for  the  poor 
Agreement . 


INTRODUCTION 

A  knowledge  of  viscous  effects  in  transonic  and  supersonic 
flows  has  become  increasingly  important  because  of  the  demand  for 
improved  accuracy  in  the  prediction  of  high-speed  flight  phenomena. 
In  particular,  viscous  effects  on  the  walls  of  supers onic^ channels 
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may  profoundly  affect  the  uniformity  of  the  flow  in  the  test 
stream,  thereby  imposing  limitations  on  the  usefulness  of  such 
channels  for  investigating  the  flow  on  models. 

The  potential  flow  in  a  supersonic  channel  may  be  readily 
determined  by  the  methods  of  references  1  and  2.  In  reference  3, 
a  method  is  proposed  for  predicting  the  turbulent  boundary-layer 
development  for  two-dimensional  flows  with  pressure  gradient.  In 
this  method,  the  velocity  profile  is  approximated  by  a  power  pro¬ 
file  and  the  particular  power  profile  and  the  value  of  skin-friction 
coefficient  are  obtained  from  empirical  laws  for  turbulent  boundary 
layers.  Because  these  laws  were  formulated  from  results  obtained 
at  low  speeds  (reference  4),  the  validity  of  their  extension  to 
boundary  layers  in  supersonic  streams  is  not  apparent. 

The  present  study,  which  was  conducted  at  the  NACA  Lewis  lab¬ 
oratory,  presents:  (1)  an  evaluation  of  the  method  of  reference  3 
for  predicting  the  boundary-layer  development  along  the  contour 
walls  and  the  side  walls  of  a  typical  supersonic  channel;  and  (2) 
an  investigation  of  the  validity  of  the  power-prof ile  parameter 
and  the  skin-friction  coefficient  assumed  in  the  method.  The 
actual  flow  in  a  two-dimensional  supersonic  channel,  however,  dif¬ 
fers  from  the  idealized  flow  assumed  in  the  analysis  in  reference  3 
because  of  comer  effects  and  secondary  flows;  hence  a  precise  check 
between  experiment  and  theory  should  not  be  expected.  An  attempt 
was  therefore  made  to  account  for  the  differences  between  experi¬ 
ment  and  theory  in  terms  of  secondary-f  low  phenomena  and  to  eval¬ 
uate  these  effects  qualitatively  in  the  channel  investigated. 


APPARATUS 

A  diagram  of  the  3.84-  by  10-inch  test  channel  used  in  this 
investigation  with  the  adjacent  surge  chamber,  diffuser,  piping, 
valves,  and  screen  and  pertinent  dimensions  is  shown  in  figure  1. 
The  maximum  pressure  attainable  at  the  inlet  was  40  inches  of  mer¬ 
cury  absolute  and  the  minimum  exhaust  pressure  was  l|*  inches  of 
mercury  absolute. 

The  supersonic  nozzle  was  designed  to  give  a  potential -flow 
Mach  number  of  2.08  by  the  method  of  reference  1  with  a  throat  of 
conservatively  large  radius  and  is  shown  to  scale  in  figure  2. 

The  bellmouth  contraction,  the  Contour  of  which  is  not  shown  in 
figure  2,  consisted  of  smoothly  faired  surfaces  generating  an  area 
contraction  ratio  of  12. 
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The  entire  channel  from  the  bellmouth  inlet  to  the  subsonic 
diffuser  was  constructed  of  stainless  steel  machined  and  polished 
to  a  maximum  surface  irregularity  of  10  microinches.  Nozzle  coor¬ 
dinates  were  accurate  to  ±0.010  inch. 

The  locations  at  which  the  boundary-layer  total -pressure  sur¬ 
veys  were  made  are  also  shown  in  figure  2.  Static -pres  sure  ori¬ 
fices  0.013  to  0.020  inch  in  diameter  were  placed  at  these  locations 
and  four  thermocouples  were  embedded  in  the  side  wall  near  the  sur¬ 
face  at  the  positions  indicated. 

Both  single-tube  and  rake  probes  were  used  for  measuring  the 
boundary -layer  profiles.  The  single -tube -probe  tips  were  made  of 
stainless-steel  tubing  having  a  0.020-inch  outside  diameter  flat¬ 
tened  to  0.006-inch  outside  thickness.  The  rake-probe  tips  varied 
from  0.015-  to  0.050-inch  diameter.  Some  of  the  0.015-inch- 
diameter  tips  were  flattened  to  0.005 -inch  outside  thickness, 
others  to  0.010,  and  still  others  were  left  circular;  the  flattened 
tips  were  nearest  the  wall. 

The  length  of  the  probes  was  determined  by  noting  how  far 
upstream  the  wall  static  pressure  was  disturbed  when  a  dummy -probe 
support  was  introduced  into  the  stream.  The  probe  was  then  designed 
to  place  the  tip  considerably  upstream  of  this  point.  The  four 
probes  used  in  this  investigation  are  shown  in  figure  3.  Probes  a 
and  b  were  3  inches  long,  c  was  4  inches,  and  d  was  variable 
from  6  to  12  inches  .  Probe  d  did  not  possess  sufficient  rigidity 
for  great  accuracy  but  the  length  was  necessary  for  transonic  and 
subsonic  measurements. 


PROCEDURE 

Boundary -layer  surveys  were  made  in  the  positions  indicated 
in  figure  2  using  probes  a,  b,  and  d  on  the  side  wall  and 
probe  c  on  the  bottom  wall.  Probe  d  was  used  at  x  <.  5.3  inches, 
prohe  b  at  x  >9.3  inches,  and  probe  a  at  several  positions 
to  check  the  measurements  of  probe  b,  where  x  is  the  distance 
in  inches  downstream  of  the  geometrical  throat  of  the  nozzle. 

For  each  x  location  of  the  prohe,  the  stagnation  pressure 
was  varied  from  37  to  5  inches  of  mercury  absolute  for  most  of  the 
runs,  although  some  of  the  results  presented  range  from  40  to 
5  inches.  The  increments  between  the  various  inlet  pressures 
investigated  varied  from  4  to  24  inches  of  mercury;  the  smaller 
increment  was  used  in  the  pressure  range  at  which  appreciable 
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changes  in  the  boundary-layer  development  took  place.  The  inlet 
temperature  and  dew-point  temperature  at  atmospheric  pressure  were 
selected  as  130°  and  -20°  F,  respectively, . to  give  condensation- 
free  flow  at  all  inlet  pressures .  Pressure  data  were  taken  only 
after  the  channel  side -wall  temperature  reached  equilibrium,  as 
indicated  by  the  embedded  thermocouples . 

The  conventional  equations  for  reducing  pressure  data  to  Mach 
number  using  subsonic  and  supersonic  pitot  tubes  were  used  (refer¬ 
ence  2,  pp.  26  and  77). 

EXPERIMENTAL  ACCURACY 

The  maximum  total  variations  of  the  stagnation  conditions  dur¬ 
ing  a  given  boundary-layer  survey  or  from  one  run  to  another  from 
the  conditions  previously  stated  were: 


Reservoir  pressure,  in.  Hg . ±0.05 

Reservoir  temperature,  °F . .  ±3 

Dew  point  at  atmospheric  pressure,  °F  ..........  .  ±8 


In  evaluating  the  accuracy  of  the  boundary -layer  total- 
pressure  measurements ,  three  sources  of  possible  error  must  be 

considered:  (1)  manometer  error,  (2)  probe -positioning  error,  and  * 

(3)  probe  influence  on  the  boundary  layer. 

As  a  conservative  value,  the  maximum  manometer  error  was 
estimated  at  a  consistent  ±0.1  inch  of  mercury  for  a  given  survey 
and  had  its  greatest  effect  on  the  computed  values  of  M,  8*, 
and  6  when  applied  to  the  static -pres  sure  measurement.  (The  sym¬ 
bols  used  herein  are  defined  in  the  appendix.)  The  chief  probe - 
positioning  error  was  in  zeroing  the  probe  against  the  channel  wall 
and  had  a  consistent  maximum  value  of  ±0.002  inch  for  probes  a, 
b,  and  c  and  ±0.005  inch  for  probe  d  for  a  given  survey.  At 
inlet  pressures  of  37,  13,  and  5  inches  of  mercury  absolute  and 
Mach  number  of  2.0,  the  maximum  errors  due  to  manometer  and  posi¬ 
tioning  errors  are  summarized  in  the  following  table: 


Inlet 

Error 

Error 

Error 

Error 

pressure 

p0 

(in.  Hg) 

in  M 

(percent) 

in  5 

( percent ) 

in  5* 
(percent) 

in  8 

(percent) 

37 

2.3 

0.5 

2.0 

13 

6.7 

.5 

2.6 

5 

17.3 

.5 

6.7 
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For  stations  upstream  of  the  test  section,  where  the  static  pres¬ 
sure  increases,  the  percentage  error  in  measuring  Mach  number 
should  he  less  than  the  values  shown  in  the  preceding  table  for  a 
given  inlet  pressure.  The  percentage  error  in  5*  and  0,  how¬ 
ever,  should  increase  as  the  probe  is  moved  upstream  because  the 
increased  positioning  error  associated  with  reduced  boundary-layer 
thicknesses  will  overbalance  decreased  errors  in  measuring  Mach 
number. 

Other  errors,  such  as  those  resulting  from  the  influence  of 
the  probe  on  the  boundary  layer,  could  not  be  checked  quantita¬ 
tively.  The  observation  was  made,  however,  that  static  pressures 
increased  up  to  2  percent  as  the  probe  tip  approached  the  wall. 
This  effect  was  present  with  probes  of  length  3  to  12  inches  and 
could  not  be  eliminated  by  lengthening  the  probe. 

A  comparison  of  boundary-layer-profile  measurements  using  the 
single-tube  probe  (probe  a)  and  the  rake  (probe  b)  showed  excel¬ 
lent  agreement  in  measured  values  of  Mq  and  5.  A  random  varia¬ 
tion  in  the  profile  shape’  from  one  test  to  another  obtained  from 
both  probes,  however,  produced  a  variation  of  0  to  10  percent  in 
the  value  of  5*.  The  effective  probe  center  was  considered  to  be 
at  the  geometrical  center  of  the  probe  tip. 

The  accuracy  of  the  temperature  measurements  on  the  side  wall 
was  ±2°  F. 


RESULTS  AMD  DISCUSSION 
Mach  Number  Distributions  along  Channel 

A  requisite  for  evaluating  the  experimental  flow  measurements 
of'  this  investigation  is  that  the  theoretical  potential -flow  Mach 
number  distribution  throughout  the  channel  be  known.  In  order  to 
obtain  the  theoretical  Mach  number  distributions  along  the  curved 
bottom-wall  and  the  flat  side-wall  center  lines  in  the  supersonic 
part  of  the  channel,  the  method  of  characteristics  (reference  2) 
using  the  assumption  of  uniform  parallel  flow  in  the  nozzle  throat 
was  employed.  These  distributions  are  indicated  in  figure  4. 

Shown  also  in  figure  4  are  the  experimental  stream  Mach  num¬ 
ber  distributions  for  a  range  of  inlet  pressures  Pq  of  5  to 
40  inches  of  mercury  absolute.  The  Mach  number  was  computed  by 
two  methods;  (1)  from  measured  local  total  and  static  pressures, 
and  (2)  from  measured  inlet  total  and  local  static  pressures.  Mach 


6 


MCA  m  2203 


numbers  in  the  test  section  computed  by  the  first  method  decreased 
on  the  bottom  -wall  (fig.  4(a))  as  the  inlet  pressure  was  reduced 
and  increased  on  the  side  wall  (fig.  4(b)).  Conversely,  the  Mach 
numbers  computed  from  measured  inlet  total  and  local  static  pres¬ 
sures  with  no  probes  in  the  stream  increased  on  both  walls  with 
decreasing  pressures.  No  reason  for  this  behavior  is  known. 

Further  comparisons  of  the  Mach  number  at  Pq  values  of  13 
and  5  inches  of  mercury  (fig.  4)  showed  large  irregularities  in 
the  distribution  when  Mach  number  was  computed  by  the  first  method 
and  smooth  distributions  when  using  the  second.  Part  of  the  irreg¬ 
ularity  in  using  the  first  method  was  no  doubt  caused  by  decreased 
accuracy  in  the  low-pressure  results,  but  the  greater  portion  of 
it  is  believed  to  be  due  to  actual  irregularities  in  the  air  flow, 
possibly  weak  shock  waves.  These  irregularities  in  the  flow  did 
not  affect  the  static  pressure  at  the  wall  and  hence  are  not  indi¬ 
cated  in  the  Mach  number  computed  by  assuming  a  constant  stagnation 
pressure . 

The  Mach  number  distributions  obtained  from  the  inlet  total 
and  local  static  pressures  were  included  in  figure  4  for  comparison 
only  and  were  not  used  in  the  analysis  of  the  results .  The  curves 
of  theoretical  potential  flow,  corrected  for  experimental  boundary- 
layer  development,  will  be  discussed  in  the  section  "Boundary- 
Layer  Development . ” 

Another  observation  to  be  made  from  figure  4  is  that  the 
experimental  Mach  numbers  obtained  from  local  total  and  static 
pressure  frequently  exceed  the  potential -flow  values  at  Pq  of  13 
and  5  inches  of  mercury.  Although  the  amounts  that  these  Mach  num¬ 
bers  exceed  the  potential -flow  values  are  very  close  to  the  maximum 
experimental  error,  the  consistency  of  these  results  and  the 
appearance  of  the  same  results  in  independent,  unpublished  tests 
corroborate  the  present  findings.  Also,  the  values  on  the  side 
wall  are  much  higher  than  those  on  the  bottom  wall  at  these  low 
pressures  for  a  large  part  of  the  channel  length. 

Lastly,  a  comparison  of  the  Mach  number  distributions  for 
Pq  of  37  and  21  inches  of  mercury  (fig.  4(b))  shows  them  to  be 
practically  the  same,  which  means  that  no  significant  changes  in 
the  potential -flow  distribution  or  boundary -layer  development 
occurred  in  this  pressure  range. 
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Temperature  Recovery  and  Prandtl  Number 


Experimental  recovery  factors  r\r  were  computed  from  measured 
reservoir  temperature  Tq,  measured  channel -vail  surface  tempera¬ 
ture  Tw>  and  computed  stream  temperature  Tq  based  on  local 
experimental  Mach  number  and  were  defined  by 


These  values  of  T]r  were  then  compared  with  theoretical  values 
for  an  insulated  plate  given  by  the  one-third  power  of  the  Prandtl 
number  Pr  (reference  5),  where  the  value  of  -Pr  was  a  function 
of  the  temperature  at  the  wall  (reference  6).  The  following  table 
summarizes  these  results: 


X 

(in.) 

Mi 

(°F) 

% 

1/3 

Pr 

-1.3 

0.88 

118 

0.852 

0.888 

18.3 

1.76 

102 

.876 

.889 

36.3 

2.04 

98 

.881 

.890 

54.3 

2.01 

98 

.879 

.890 

At  the  first  position,  where  the  wall  temperature  was  about  40°  F 
above  room  temperature,  the  value  of  qr  was  about  4  percent 
l/3 

below  Pr  '  ;  whereas  at  the  three  remaining  positions,  for  which 
the  wall  temperature  was  only  about  20°  F  above  room  temperature, 

the  value  of  T)r  was  less  than  ljy  percent  below  Pr1^3 .  These 

results  indicate  that  heat  transfer  through  the  walls  is  small  and 
may  probably  be  neglected  in  computing  velocities  in  the  boundary 
layer  and  other  boundary -layer  quantities. 

The  assumption  of  zero  heat  transfer  made  in  reference  3  and 
used  in  the  present  analysis  has  been  shown  to  be  approximately 
correct,  but  a  question  still  remains  concerning  the  error  intro¬ 
duced  by  considering  the  Prandtl  number  equal  to  1  in  the  interpre¬ 
tation  of  the  data.  This  problem  is  analyzed  in  detail  in  refer¬ 
ence  7,  where  it  was  found  that  the  use  of  Pr  =  1  led  to  an 
error  of  3  percent  in  the  velocity  near  the  wall  surface,  which 
rapidly  diminished  to  0  as  the  free  stream  (M  =  2.0)  was  approached. 
As  a  result  of  these  errors  in  velocity,  subsequent  errors  of 
0.1  percent  were  introduced  in  5*  and  9.  The  errors  introduced 
into  the  skin-friction  measurements  also  were  small.  The  simplifying 
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assumptions  of  zero  heat  transfer  and  a  Prandtl  number  of  1  used 
in  the  following  analysis  therefore  appear  to  be  justified. 


Boundary -Layer  Development 

Experimental  boundary -layer  development.  -  Nondimeris ional 
boundary-layer  velocity  profiles  u/u$  against  j/E  for  the  side 
wall  in  the  transonic -flow  region  of  the  channel  (nozzle  throat) 
at  various  inlet  pressures  are  plotted  in  figure  5.  The  boundary- 
layer  thickness  5  is  defined  as  that  distance  in  the  y -direction 
(perpendicular  to  the  wall)  at  which  the  velocity  u §  =  0.99  u^, 
where  U]_  is  the  asymptotic  stream  velocity.  These  profiles  were 
measured  throughout  a  range  of  free -stream  Mach  numbers  from  0.57 
to  1.68.  The  profiles  at  Pq  of  37.  and  13  inches  of  mercury  are 
typical  turbulent  profiles  throughout;  whereas  those  at  Pq  of 
5  inches  are  laminar  for  x  <  5.3  inches  and  turbulent  for 
x  >  9.3  inches.  Theoretical  laminar  profiles  (reference  8)  assum¬ 
ing  no  heat  transfer  are  also  presented  for  the  lowest  inlet  pres¬ 
sure  corresponding  to  the  experimental  Mach  numbers  at  x  of  -6.7 
and  5.3  inches.  Because  the  theoretical  laminar  profiles  were  based 
on. flat-plate  theory,  whereas  the  experimental  layers  appeared  in 
a  highly  favorable  pressure  gradient,  the  disagreement  between  the 
two  should  not  be  surprising.  The  similarity  of  the  slopes  near 
the  wall,  however,  identifies  the  experimental  profiles  as  laminar. 

Nondimens ional  velocity  profiles  of  the  boundary  layer  through¬ 
out  the  entire  channel  are  presented  in  figure  6  at  Pq  of  37  inches 
of  mercury  for  both  the  bottom  and  side  walls.  Logarithmic  coor¬ 
dinates  were  used  to  facilitate  comparisons  between  the  theoretical 
power-law  profiles  and  the  experimental  points.  The  theoretical 
profiles  are  represented  by  the  straight  lines,  with  the  appropriate 
value  of  the  power-profile  parameter  N,  which  defines  the  shape 
of  the  theoretical  boundary -layer  power  profile  given  by 


In  comparing  the  experimental  and  theoretical  profiles,  it  is  seen 
that  the  best  agreement  is  reached  at  the  greater  distances  down¬ 
stream,  which  correspond  to  the  regions  of  near  zero  pressure  gra¬ 
dient.  The  results  shown  in  figure  6  will  be  further  analyzed  in 
a  following  discussion  of  the  power-profile  parameter. 

The  boundary -layer  development  in  terms  of  the  boundary -layer 
thickness  5,  the  displacement  thickness  5*,  and  the  momentum 
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thickness  0  for  the  bottom  and  side  walls  at  inlet  pressures 
from  37  to  5  inches  of  mercury  absolute  are  presented  in  figure  7. 
Also  included  are  curves  faired  near  the  experimental  points  that 
were  used  in  evaluating  skin-friction  coefficients. 

The  boundary-layer  developments  for  Pq  of  37  inches  of  mer¬ 
cury  on  the  bottom  and  side  walls  are  characterized  by  smooth  dis¬ 
tributions  of  5,  5*,  and  6  along  the  channel  length.  The 

development  for  Pq  of  21  inches  of  mercury  along  the  side  wall 
(not  shown  in  fig.  7(b))  was  nearly  the  same  as  that  obtained  for 
Pq  of  37  inches  of  mercury.  Decreases  in  inlet  pressure  below 
21  inches  of  mercury  resulted  in  increases  in  the  rate  of  boundary- 
layer  growth  in  the  turbulent  regions  and  decreases  in  the  laminar 
regions.  In  addition,  the  boundary -layer  developments  along  the 
bottom  wall  and  the  Mach  number  distributions  of  figure  4(a)  (com¬ 
puted  from  local  total  and  static  pressures)  become  highly  irreg¬ 
ular.  These  irregularities  in  the  boundary -layer  development  are 
not  caused  primarily  by  the  decreased  accuracy  of  the  low-pressure 
measurements,  but  rather  reflect  the  actual  irregularity  of  the 
boundary -layer  thickness,  the  measurement  of  which  does  not 
strongly  depend  on  the  accuracy  of  the  pressure  measurements. 

Another  characteristic  of  the  low  inlet -pressure  results  is 
the  apparent  incompatibility  of  the  Mach  number  distributions  along 
the  channel  with  the  boundary-layer  development.  Boundary -layer 
displacement  thicknesses  are  considerably  greater  for  the  low- 
pressure  flow  downstream  in  the  channel  and  less  in  the  throat 
region;  yet  substantial  increases  in  Mach  number  were  noticed  as 
the  inlet  pressure  was  reduced.  This  anomaly  was  previously 
pointed  out  when  it  was  shown  that  the  experimental  Mach  number 
exceeded  the  theoretical  potential -flow  Mach  number  at  low  pres¬ 
sures.  Comparison  of  the  theoretical  Mach  number  distribution, 
which  has  been  corrected  for  the  presence  of  the  experimental  dis¬ 
placement  thickness,  with  the  experimental  Mach  number  distribution 
throughout  the  range  of  inlet  pressures  emphasizes  this  anomaly. 
Curves  of  theoretical  Mach  number  corrected  for  6*  are  given  in 
figure  4  for  the  bottom  and  side  walls.  In  correcting  the  potential- 
flow  Mach  number  it  was  assumed  that  transverse  to  the  stream 
direction  the  boundary  layer  was  of  constant  thickness  and  veloc¬ 
ity  distribution.  The  geometric  channel-area  ratios  were  then  cor¬ 
rected  for  the  displacement  thickness  and,  with  the  assumption  of 
one -dimensional  flow,  the  corrected  potential  flow  Mach  numbers 
were  obtained . 

A  comparison  of  the  corrected  potential  flow  and  experimental 
Mach  number  distributions  in  figure  4  indicates  better  agreement 
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for  the  high-inlet-pressure  range  than  for  the  low  pressures.  As 
Po  is  progressively  reduced  to  13  and  5  inches  of  mercury ,  increas¬ 
ing  disagreement  is  noted,  that  is,  the  corrected  potential -flow 
Mach  number  is  generally  less  than  the  experimental.  The  boundary- 
layer  displacement  thickness  is  believed  to  vary  in  a  direction 
transverse  to  the  stream,  thereby  invalidating  the  one -dimensional 
area-ratio  correction.  The  presence  of  such  nonuniformity  of  the 
boundary  layer  is  possible  if  secondary  flows  are  present.  The 
variation  in  the  discrepancy  between  the  experimental  and  corrected 
potential -flow  Mach  numbers  in  figure  4  as  the  inlet  pressure  is 
changed  suggests  a  variation  in  the  intensity  of  the  secondary 
flows  with  Reynolds  number. 


Reynolds  number  and  transition.  -  In  reference  3,  two  defini¬ 
tions  of  Reynolds  number  were  investigated  in  the  determination  of 
the  turbulent  skin-friction  coefficient  in  a  pressure  gradient. 

One  definition  is  based  on  the  kinematic  viscosity  of  the  free 
stream  and  the  other,  on  the  kinematic  viscosity  of  the  air 

adjacent  to  the  wall  in  the  boundary  layer  uw.  They  may  be  writ¬ 
ten  as 


Re—  n  = 

X,1 


ulxl 


(3a) 


and 


Re; 


x.w 


unx 
1  w 

V 

w 


(3b) 


The  terms  x^  and  5c^  are  the  equivalent  lengths  of  turbulent 
run  on  a  flat  plate  necessary  to  produce  a  known  boundary -layer 
momentum  thickness.  These  lengths  of  run  were  computed  assuming  a 
flat-plate  skin-friction  coefficient  dependent  on  Reynolds  number 
based  on  stream  and  wall  kinematic  viscosities,  respectively. 


The  preceding  methods  for  evaluating  Reynolds  number  are 
indirect  and  may  be  replaced  by  more  convenient  definitions  when 
6  is  known;  such  as 


Re 


(4a) 


and 


R9e,w 


Uf  6 


(4b) 


t 


r 
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The  applicability  of  these  two  definitions  of  Reynolds  number  will 
*  be  investigated  further  in  the  discussion  of  the  shin-friction 

coefficient. 

Laminar  boundary  layers  were  observed  at  Pq  of  5  inches  of 
mercury  for  values  of  x  5.3  inches  on  the  side  wall.  At 
x  of  9.3  inches  (the  next  position  investigated)  and  at  greater 
distances  downstream,  turbulent  layers  were  found.  Reynolds  num¬ 
bers  Reg  3_  and  Reg  w  and  those  based  on  5*,  Re^*  p  and 
Re§*  v  at  x  of  5.3  and  9.3  inches  are  listed  with  values  given 
by  Goldstein  (reference  9)  for  the  commencement  of  transition  in 
the  following  table : 


Distance  from 
throat,  x 
(in.) 

Transition 
(Goldstein, 
reference  9) 

5.3 

Laminar 

9.3 

Turbulent 

Ee0,l 

288 

845 

Ee0,V 

193 

433 

Ee5*,l 

913 

1982 

560  to  1700 

Ee6*,w 

612 

1015 

Apparently  transition  took  place  near  the  range  given  by  Goldstein 
and  either  Re§*  ^  or  Re§*  v  falls  into  it.  The  lowest  value  of 
Reg*  i  computed  for  the  bottom  wall  was  1160  and  no  laminar 
boundary  layers  were  observed  there. 

Form  parameters.  -  In  reference  10,  von  Doenhoff  and  Tetervin 
suggest  the  usefulness  of  the  turbulent -form  parameter  Hj_  for 
defining  the  shape  of  the  boundary-layer  velocity  profile  at  low 
speeds  where  compressibility  is  negligible.  In  the  parameter 
Hj_  =  5^*  and  6^  are  the  displacement  and  momentum 

thicknesses,  respectively,  which  are  computed  assuming  constant 
stream  density  in  the  boundary  layer  normal  to  the  wall. 

Turbulent -form  parameters  computed  assuming  a  variable  density 
in  the  boundary  layer  H  and  form  parameters  are  presented 

in  figure  8  for  the  bottom  and  side  walls  at  the  various  inlet 
pressures .  An  outstanding  characteristic  of  these  plots  is  the 
constancy  of  the  values  of  Hi  of  1.20  for  the  bottom  and  1.24 
for  the  side  wall  throughout  the  supersonic  region,  with  the 


12 


MCA  TO  2203 


exception  of  the  throat  region,  which  suggests  that  the  velocity 
profile  may  be  independent  of  the  Mach  number.  The  increases  in 
H  and  Hi  near  the  throat  were  the  result  of  changes  in  velocity 
profile;  the  irregularities  in  the  distributions  there  were  caused 
by  a  greater  uncertainty  inherent  in  the  boundary-layer  measure¬ 
ments  because  of  the  very  small  boundary-layer  thicknesses  and 
large  probe  flexibility. 

Also  included  in  figure  8  are  theoretical  values  of  H 
obtained  using  the  theory  of  reference  3.  The  agreement  between 
theory  and  experiment  appears  to  be  better  at  the  high  inlet  pres¬ 
sures  where  the  boundary -layer  development  was  smooth  than  at  the 
low  pressures,  and  appears  better  at  values  of  x  greater  than 
8  inches  than  at  values  less  than  8  inches.  Values  of  H  ranged 
from  about  1.5  near  the  channel  entrance  to  about  3.2  farther  down¬ 
stream  for  both  experiment  and  theory.  This  large  range  of  values 
is  directly  attributable  to  the  compressibility  of  the  air  at  high 
Mach  numbers;  the  large  magnitude  of  H  at  the  downstream  posi¬ 
tions  is  attributable  to  the  low  density  of  the  air  in  the  boundary 
layer  near  the  wall.  Hence,  the  predominant  change  in  the  boundary - 
layer-profile  development  occurred  in  the  density  distribution. 

Comparison  of  experimental  and  theoretical  boundary -layer 
development .  -  A  theoretical  method  for  determining  the  turbulent 
boundary -layer  development  using  the  Karman  momentum  equation  is 
given  in  reference  3  and  is  used  to  predict  the  boundary -layer 
development  in  the  two-dimensional  supersonic  channel  investigated 
in  this  report.  The  following  assumptions  in  the  method  are 
discussed: 

1.  No  heat  is  transferred  through  the  channel  walls  and  use 
of  a  Prandtl  number  equal  to  1  does  not  lead  to  inadmissible  error 
in  computing  the  boundary-layer  development. 

2.  Empirical  laws  for  the  turbulent  boundary-layer  power- 
profile  parameter  and  skin-friction  coefficient  obtained  at  low 
subsonic  speeds  in  zero  pressure  gradients  can  be  extended  to 
supersonic  flows  in  highly  favorable  pressure  gradients. 

3.  Theoretical  boundary-layer  development  does  not  affect  the 
theoretical  potential  flow  sufficiently  to  make  a  second  calcula¬ 
tion  of  the  boundary -layer  development  necessary. 

4.  The  pressure  gradients  transverse  to  the  stream  direction 
do  not  produce  appreciable  secondary -flow  effects . 
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The  first  of  these  assumptions  has  already  been  considered 
and  has  been  shown  to  have  a  negligible  effect  on  the  boundary- 
layer  development  and  the  shin-friction  coefficient.  The  accuracy 
in  the  remaining  assumptions  taken  collectively  may  be  estimated 
by  a  comparison  of  the  experimental  and  theoretical  boundary-layer 
developments . 

In  f  igure  9,  the  experimental  turbulent  boundary-layer  develop¬ 
ments  in  terras  of  8,  5*,  and  0  at  Pq  of  37  inches  of  mercury 
as  obtained  from  experiment  and  theory  (reference  3)  are  presented. 
Two  theoretical  developments  were  computed  using  turbulent  power- 
profile  parameters  and  skin-friction-coefficient  laws  based  on 
He£  i  and  Be“  v.  A  comparison  of  these  two  theoretical  develop¬ 
ments  with  the  experimental  development  shows  that  the  theoretical 
development  in  which  Re—  v  is  used  to  estimate  the  power-profile 
parameters  and  skin-friction  coefficients  is  the  better  by  far  on 
both  the  bottom  and  side  walls.  In  almost  every  case,  the  use  of 
Ee~  -j_  gives  theoretical  boundary  layers  that  are  far  thicker  than 
experimental  throughout  the  channel. 

Because  of  it£  importance  in  determining  the  skin  friction 
from  profile  measurements,  the  rate  of  growth  of  the  boundary 
layer  expressed  by  8,  8*,  and  6  must  be  considered.  If  the  rate 
is  included  in  the  comparison  it  will  be  noted  that  even  in  the 
case  of  the  theoretical  boundary -layer  development  computed  using 
Be-  v  significant  disagreements  occur  between  experiment  and  the¬ 
ory.  This  trend  is  especially  true  on  the  side  wall  where  the 
experimental  growth  near  the  entrance  to  the  channel  (x  <:  20  in.) 
is  less  than  theoretical  and  farther  downstream  where  it  is 
greater.  On  the  bottom  wall,  the  reversed  trend  is  apparent  but 
to  a  lesser  degree.  These  differences  in  the  rate  of  boundary- 
layer  growth  are  probably  caused  by  secondary  flows  in  the  boundary 
layer  and  will  be  discussed  more  adequately  in  a  later  section  on 
skin-friction  coefficient. 

A  second  calculation  of  the  boundary -layer  development  was 
made  to  determine  the  effect  of  the  theoretical  displacement 
thickness  on  the  assumed  potential  flow.  A  one -dimensional  cor¬ 
rection  similar  to  the  one  employed  earlier  in  correcting  the 
potential -flow  Mach  number  distribution  for  the  presence  of  the 
experimental  displacement  thickness  was  used;  the  resulting  cor¬ 
rected  potential  flow  was  used  to  recompute  the  boundary -layer 
development  a  second  time.  Changes  in  the  recomputed  development 
were  less  than  2  percent,  indicating  sufficient  accuracy  in  using 
the  first  approximation  only. 


14 


NACA  TN  2203 


Power-profile  parameter.  -  The  power-prof ile  parameter  E, 
■which  defines  the  shape  of  the  boundary -layer  profile  according  to 
equation  (2),  was  evaluated  from  the  experimental  Mach  number  and 
turbulent -form-parameter  distributions  along  the  channel  length 
for  an  inlet  pressure  of  37  inches  of  mercury.  The  distribution 
of  N  thus  obtained  is  presented  in  figure  10  for  the  bottom  and 
side  walls.  In  order  to  obtain  E  from  experimental  values  of  M  ^ 
and  H,  use  was  made  of  table  III  of  reference  3,  in  which  N  is 
given  as  a  function  of  Mj_  and  E.  Because  N  is  extremely  sensi¬ 
tive' to  small  variations  in  either  Mpor  H  and  because  the  var¬ 
iation  of  H  with  x  in  figure  8  is  not  smooth,  faired  values  of 
H  were  used  in  determining  E. 


Empirical  values  of  the  power-profile  parameter  are  also  pre¬ 
sented  in  figure  10  to  be  compared  with  the  experimental  points. 
These  empirical  values  are  given  by  the  following  formulas: 


E 


1  or 


w  "  2-6  or 


l/l4 


which  was  used  in  reference  3,  by  setting  u/u^  =  1.  Very 
difference  between  Ej  and  is  apparent  for  the  range 
Reynolds  number  encountered  in  this  experiment  because  the 
nent  l/l4  is  small;  hence  either  Ep  or  %  is  equally 
able  here. 


little 

of 

expo- 

accept- 


A  comparison  between  the  experimental  and  empirical  power- 
profile  parameters  shows,  in  general,  that  the  empirical  values 
underestimate  the  experimental  values  on  the  bottom  wall  and  over¬ 
estimate  them  on  the  side  wall.  In  consideration  of  only  the 
proximity  of  the  side  walls,  it  would  appear  that  the  experimental 
power-profile  parameter  along  the  bottom-wall  center  line  would  be 
less  than  along  the  side-wall  center  line  because  the  boundary  layer 
along  the  bottom  wall  is  subject  to  a  large  extent  to  the  viscous 
action  produced  by  the  side  walls  near  the  channel  corners.  The 
actual  behavior  of  the  power-profile  parameter,  however,  contradicts 
this  supposition;  the  only  explanation  for  the  actual  behavior  appears 
to  be  in  the  existence  of  secondary  flows . 

Another  method  for  determining  N  for  the  experimental  pro¬ 
files  would  be  to  find  directly  the  value  of  E  that  corresponds 
most  closely  to  the  experimental  profiles  in  figure  6.  A  selection 
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of  a  value  of  N  corresponding  to  the  experimental  profile  is  dif¬ 
ficult,  however,  because  the  experimental  points  actually  cross 
several  values  of  N,  particularly  for  the  upstream  boundary 
layers.  Only  for  the  downstream  profiles  can  a  rational  selection 
of  N  be  made,  and  there  the  agreement  between  this  new  value  of 
I  with  the  previously  determined  experimental  value  will  be  good. 


Skin-Friction  Coefficient 


Experimental  skin-friction  coefficients  (cf)q  were  computed 
along  the  wall  center  lines  by  substitution  of  the  faired  values 
of  M]_,  6*,  and  6  presented  in  figures  4  and  7  into  the  Karman 
momentum  equation  (reference  3)  for  two-dimensional  flow  along  the 
bottom  wall 


d6 

dx 


+ 


6  (2-M-j2)  +  5* 

Mi  ( 1  ♦  ¥  «i2) 


dMx 

dx 


fw  T 

Pi  Pwul2 


(°f  )l 


(6a) 


and  three-dimensional  divergent  flow  along  the  side  wall 


de 

6  (2-Mq2)  +  5* 

dr  + 

«i  ( 1  -  ¥  Wi2) 

dM,  6 
dr  r  = 


'w 


P1  pwUl 


(cf)1  (6b) 


The  experimental  friction  coefficients  so  computed  are  presented 
in  figures  11  to  13  for  the  bottom  and  side  walls  at  Pq  of  37 
to  5  inches  of  mercury. 

Accurate  calculation  of  skin-friction  coefficients  by  the  pre¬ 
ceding  method  is  difficult  particularly  at  the  low  inlet  pressures 
where  considerable  judgment  is  necessary  for  fairing  curves  for 
the  Mach  number  distribution  and  boundary -layer  development.  The 
measurement  of  derivatives  in  equations  (6a)  and  (6b)  makes  this 
task  even  more  difficult.  For  these  reasons,  it  was  necessary  to 
use  curves  that  were  faired  near  the  experimental  values  of  M 
and  the  boundary -layer  development  and  to  avoid  inflection  points 
except  in  the  transition  region  of  the  boundary  layer  on  the  side 
wall  at  Pq  of  5  inches  of  mercury .  Approximate  average  devia¬ 
tions  in  (CfO-L  computed  from  faired  and  unfaired  results  are 
listed  in  the  following  table.- 
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Inlet 

pressure 

Deviations  in  (cf)i 
( percent ) 

P0 

Bottom 

Side 

(in.) 

vail 

vail 

37 

±3 

±10 

13 

±10 

±20 

5 

±20 

±20 

Friction  coefficients  along  the  'bottom  vail  dropped  severely 
from  values  averaging  0.0023  at  x  =  0  for  the  three  inlet  pres¬ 
sures.  In  contrast,  more  gradual  drops  beginning  at  (cf)i  of  0.0011 
to  0.0016  took  place  on  the  side  vails  at  x  =  0.  Generally,  fric¬ 
tion  coefficients  on  the  bottom  vail  continued  to  drop  smoothly  to 
x  =  44  inches  (the  last  position  surveyed  on  the  bottom  vail); 
the  variation  along  the  side  vail  vaa  more  erratic.  On  the  side 
vail  a  pronounced  increase  in  (cf).i  occurred  beginning  at 
x  =  16  inches  continuing  to  x  =  64  inches  for  Pq  of  37  inches 
of  mercury.  At  P0  of  13  inches,  the  side-vall  skin  friction 
tends  to  remain  fairly  constant  initially  vith  a  final  sving  upward 
farther  dovnstream.  The  variation  of  (cf)i  ^0  oP  ^  inches 

of  mercury  shows  the  peak  in  (cf)i  near  x  =  5  inches  taking  place 
in  the  transition  region  and  leveling  off  similarly  to  the  distri¬ 
bution  for  Po  of  13  inches  of  mercury,  but  somewhat  higher. 

Apparent  breaks  in  the  distribution  of  (cf)^  along  the  side 
vails  at  all  three  inlet  pressures  near  x  =  32  inches  are  directly 
attributable  to  the  fact  that  the  boundary-layer  flow  does  not  fol¬ 
low  the  assumed  potential -flow  streamlines.  Such  abrupt  changes 
in  (cf)x  obviously  cannot  exist  and  arise  simply  from  elimination 
of  the  term  e/r  in  equation  (6b)  at  the  vertex  of  the  test- 
section  rhombus  where  r,  as  obtained  from  the  theoretical  potential 
flow,  becomes  infinite.  Hence,  a  definite  need  exists  for  flcw- 
direction  measurements  in  boundary  layers  when  the  external  poten¬ 
tial  flow  is  rapidly  turned. 

An  empirical  turbulent  skin-friction  coefficient  lav  given 
first  by  Falkner  (reference  4)  is  restated  and  its  extension  to 
high  Mach  numbers  proposed  in  reference  3,  where  it  is  given  as 

(  \  _  0.0151 

^Cf%l  or  w 


1/7 


(7) 
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Both  (cf)x?l  (°f)x,w  were  uaeci-  to  compute  the  theoretical 

"boundary -layer  developments  presented  earlier  in  figure  9,  Falkner 
gives  an  alternate  form  of  the  preceding  law  that  can  also  be 
extended  to  depend  on  the  kinematic  viscosity  at  the  wall  as  well 
as  the  stream: 

0.00653 

=jT76  (8) 

6 , 1  or  wj 


(of  )e,i 


or  w 


When  only  stream  values  of  the  kinematic  viscosity  are  considered, 
the  respective  coefficients  are  equal: 


The  coefficients  based  on  the  kinematic  viscosity  evaluated  at  the 
wall  are  related  to  each  other  by  integrating  equation  (6a),  assum¬ 
ing  zero  pressure  gradient,  and  substituting  the  resulting  value 
of  x  in  the  expression  of  (cf)~  w,  equation  (7).  The  following 
relation  between  (cf)x,w  (°f^6,w  then  exists : 

(cf  =  (cf 


Thus  by  extending  Falkner's  formula,  which  is  based  on  U-j_,  to  for¬ 
mulas  based  on  V  three  different  values  of  the  friction  coef¬ 
ficient  can  be  defined.  These  three  values  are  presented  in  fig¬ 
ures  11  to  13  and  were  computed  from  the  experimental  Reynolds 
numbers.  In  figure  13(b),  laminar  coefficients  computed  from  ref¬ 
erence  8  are  presented  for  the  region  in  which  laminar  boundary 
layers  were  observed. 


In  presenting  the  coefficients  (cf)x  w  and  (of)e  w,  a  mul¬ 
tiplying  factor  p^/p^  was  included.  The  reason  becomes  evident 
when  the  right-hand  side  of  equations  (6a)  and  (6b)  is  shown  in 
the  following  forms: 


^w 

Pl  Pwul2 


PW  /  X 

5T  <cf>w  - 


cmi2 


(°f  )l 


The  coefficients  (cf)w  (°f)l  are  the  shin-friction  coef¬ 

ficients  evaluated  with  pw  and  p^,  respectively. 
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A  previous  comparison  in  figure  9  involving  the  experimental 
and  theoretical  boundary -layer  developments  showed  that  the  coef¬ 
ficient  (cf)“  v  gave  better  results  in  computing  the  turbulent 
boundary-layer  development  than  (cf)x,l*  A  comparison  of  (cf)x,w 
and  (cf)e  v  shows  only  slight  differences,  and  it  may  therefore 
be  concluded  that  either  of  these  coefficients  is  suitable  for 
computing  the  theoretical  turbulent  development . 

Comparisons  of  the  experimental  and  empirical  skin-friction 
coefficients  show  agreement  in  orders  of  magnitude  only  and  dis¬ 
play  considerable  differences  in  detailed  variations  along  the 
channel.  The  extremely  high  values  of  near  x  =  0  in  fig¬ 

ures  11(a)  and  12(a)  for  the  bottom  wall  are  in  sharp  contrast  to 
the  empirically  predicted  values  and  to  the  experimental  and 
empirical  coefficients  found  in  the  same  region  on  the  side  walls 
(figs.  11(b)  and  12(b)).  Likewise  the  fair  agreement  between 
experiment  and  empiricism  for  the  bottom  wall  farther  downstream 
contrasts  with  the  poor  agreement  downstream  on  the  side  wall 
where  the  experimental  coefficients  begin  to  rise  despite  the  con¬ 
tinued  increases  in  Reynolds  number. 

The  good  agreement  between  experiment  and  theory  in  the 
boundary -layer  development  using  the  empirical  skin-friction  law 
based  on  the  wall  kinematic  viscosity  (fig.  9)  as  distinguished 
from  the  poor  agreement  between  the  experimental  and  empirical 
skin-friction  coefficients  can  be  explained  on  the  basis  of  dif¬ 
ferences  in  the  rate  of  boundary -layer  growth.  One  of  the  factors 
used  in  evaluating  the  experimental  friction  coefficients  is  the 
rate  of  growth  d©/dx  and  de/dr  in  equations  (6a)  and  (6b). 

Also,  an  analysis  of  figures  9  and  11  shows  that,  for  regions  in 
which  the  experimental  boundary-layer  growth  is  more  rapid  than 
that  predicted  by  the  method  of  reference  3,  the  friction  coef¬ 
ficients  are  likewise  higher  than  empirical  and  vice  versa.  This 
analysis  considers,  of  course,  only  the  theoretical  developments 
and  friction  coefficients  involving  (°f )j  w  because  those  based 
on  stream  values  of  the  kinematic  viscosity  departed  too  far  from 
experiment.  Because  the  discrepancies  in  the  experimental  and 
theoretical  boundary-layer  growths  were  believed  to  be  the  result 
of  secondary -flow  effects  and  because  skin-friction  and  boundary- 
layer  development  are  directly  related,  the  same  explanation  would 
follow  for  the  discrepancies  between  the  experimental  and  empirical 
friction  coefficients . 
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The  theoretical  laminar  friction  coefficients  shown  in  fig- 
»  ure  13(h)  compare  reasonably  well  with  the  measured  values  in  view 

of  the  low-accuracy  characteristic  of  the  measurements  of  the  thin 
laminar  boundary  layers . 

H 

04 

CD 

Secondary  Flows 

In  both  the  theoretical  analysis  of  the  boundary-layer  devel¬ 
opment  (reference  3)  and  in  the  present  experimental  analysis,  the 
implicit  assumption  has  been  made  that  in  the  boundary  layer  no 
flow  occurs  across  planes  normal  to  the  wall  and  oriented  in  the 
direction  of  the  potential  flow,  that  is,  that  secondary  flows  are 
negligible.  This  assumption  is  known  to  be  invalid  for  flows  hav¬ 
ing  static-pressure  gradients  transverse  to  these  planes  (refer¬ 
ence  9);  the  error  introduced  by  neglecting  secondary  flows  may 
account  for  the  disagreement  between  experiment  and  theory  and  the 
various  other  anomalies  observed  in  the  analysis  of  the  experimental 
results . 

Any  nozzle  flow  will  have  the  following  characteristics: 

(a)  The  streamlines  of  the  flow  outside  the  boundary  layer  will 
tend  to  follow  along  the  predicted  potential  flow  streamlines;  and 

(b)  the  streamlines  of  the  flow  in  the  boundary -layer  will  tend  to 

«  follow  the  static -pressure  gradient.  In  all  such  nozzles,  with 

the  exception  of  axially  symmetric  ones,  these  characteristics 
wi.ll  give  rise  to  secondary  flows  in  the  boundary  layer  corre¬ 
sponding  approximately  to  the  secondary -flow  pattern  sketched  in 
figure  14.  The  static-pressure  distributions  p/pQ  along  the 
side  and  bottom  walls  causing  these  secondary  flows  are  presented 
in  the  same  figure  and  are  taken  from  the  theoretically  computed 
potential  flow. 

In  the  transonic  and  initial  expansion  region  of  the  super¬ 
sonic  nozzle,  the  boundary  layer  tends  to  flow  toward  the  top  and 
bottom  walls ;  and  in  the  straightening  region,  toward  the  center 
j.ine  of  the  side  walls.  Such  transverse  flow  of  low-energy  air 
may  account  for  some  of  the  disagreements  between  experiment  and 
theory  in  the  boundary-layer  developments  presented  earlier.  This 
flow  would  cause  the  retarded  rate  of  boundary-layer  growth  near 
the  enorance  and  the  accelerated  growth  farther  downstream  on  the 
channel -side -wall  center  line  and  the  reversed  behavior  on  the  bot¬ 
tom  wall  previously  noted.  Also,  the  low  values  of  power-profile 
parameter  N  observed  on  the  greater  part  of  the  side  wall  com¬ 
pared  with  the  higher  values  on  the  bottom  wall  suggest  a  transport 
of  low-energy  air  toward  the  center  line  of  the  side  wall. 
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More  detailed  analyses  of  other  peculiarities  observed,  such 
as  the  great  differences  in  Mach  number  between  the  side  and  bot¬ 
tom  walls  in  the  test  section  at  low  pressures,  the  presence  of 
laminar  boundary  layers  on  the  side  walls  and  turbulent  boundary 
layers  on  the  bottom  wall  near  the  channel  entrance,  and  higher 
Mach  numbers  in  the  presence  of  greater  boundary-layer  thicknesses, 
cannot  be  expected  without  specific  research  into  secondary  flows. 

The  existence  of  secondary  flows  as  previously  described  indi¬ 
cates  that  flow-direction  measurements  are  necessary  to  determine 
the  extent  of  the  three-dimensional  character  of  the  boundary -layer 
flow.  If  the  flow  departs  appreciably  from  the  assumed  two- 
dimensional  form,  momentum  equations  (6a)  and  (6b)  are  no  longer 
applicable  for  computing  skin-friction  coefficients.  For  channels 
having  higher  Mach  numbers  than  the  one  investigated  herein,  the 
•pressure  gradients  will  be  larger  and  the  secondary “flow  phenomena 
will  become  more  serious . 


.  SUMMARY  OF  RESULTS 

Boundary -layer  measurements  were  made  in  a  two-dimensional 
channel  designed  by  potential -flow  methods  for  a  uniform,  test- 
section  Mach  number  of  2.08.  Investigations  at  inlet  pressures  of  5, 
13,  and  37  inches  of  mercury  absolute  gave  the  following  results; 

1.  The  Mach  number  distribution  throughout  the  channel  and  the 
boundary -layer-prof ile  development  were  almost  independent  of  the 
inlet  pressure  in  the  range  from  37  to  21  inches  of  mercury.  From 
13  to  5  inches  of  mercury,  substantial  changes  in  the  Mach  number 
distribution  and  boundary-layer  development  took  place. 

2.  Wall-surface  temperature -recovery  factors  were  approxi¬ 
mately  1-  percent  lower  than  theoretical  recovery  factors  for 

2 

insulated  surfaces  with  turbulent  boundary  layers. 

3.  Turbulent-boundary-layer  profiles  were  found  throughout 
the  channel  at  all  inlet  pressures  with  the  exception  of  the 
entrance  region  of  the  side  wall  at  an  inlet  total  pressure  of 
5  inches  of  mercury. 

4.  The  best  approximation  of  the  boundary -layer  velocity  pro¬ 
files  to  power-law  profiles  occurred  in  regions  of  zero  pressure 
gradient . 
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5.  Transition  of  the  "boundary  layer. from  laminar  to  turbulent 
on  the  side  vail  occurred  within  the  range  of  Reynolds  number  of 
913  to  1982  based  on  boundary -layer  displacement  thickness  and 
kinematic  viscosity  evaluated  in  the  stream.  This  transition  is 
within  the  range  commonly  accepted.  No  laminar  layers  were  observed 
on  the  bottom  wall  and  Reynolds  numbers  there  were  higher  than  on 
the  side  wall. 

6.  Turbulent-f orm  parameters  computed  assuming  constant  stream 
density  in  the  boundary  layer  remained  constant  at  1.20  for  the 
bottom  and  1.24  for  the  side  wall  throughout  the  channel  for  all 
inlet  pressures  except  in  the  region  of  the  throat  where  these  val¬ 
ues  increased.  Turbulent-f orm  parameters  computed  with  variable 
density  in  the  boundary  layer  varied  from  approximately  1.5  near 
the  channel  entrance  to  3.2  in  the  test  section.  The  predominant 
change  in  the  boundary-layer-prof ile  development  occurred  in  the 
density  distribution  with  relatively  small  changes  in  the  velocity 
distribution. 

7 .  The  magnitude  of  the  experimental  boundary-layer  develop¬ 
ment  at  an  inlet  total  pressure  of  37  inches  of  mercury  agreed  well 
with  the  theoretical  magnitude  when  turbulent  skin-friction  coef¬ 
ficients  based  on  the  kinematic  viscosity  of  the  air  at  the  wall 
were  used.  The  principal  discrepancy  between  experiment  and  the¬ 
ory  was  in  the  rate  of  the  boundary-layer  development,  which  was 
probably  caused  by  secondary  flows  in  the  boundary  layer.  Because 
the  present  theory  does  not  take  secondary  flows  into  account,  a 
precise  check  between  theory  and  experiment  should  not  be  expected, 

8.  In  general,  the  experimental  power-profile  parameters  were 
higher  than  empirical  on  the  bottom  wall  and  lower  on  the  side 
wall. 

9.  Empirical  coefficients  using  Reynolds  numbers  based  on  the 
kinematic  viscosity  at  the  wall  appeared  to  conform  more  closely 
to  experiment  than  those  using  Reynolds  numbers  based  on  the  kine¬ 
matic  viscosity  in  the  stream.  A  discontinuity  in  the  distribution 
of  these  coefficients  along  the  side  wall  was  attributed  to  large 
deviations  in  the  radius  of  streamline  curvature  from  the  assumed 
potential  flow.  Such  deviations  indicate  the  need  for  flow -direction 
measurements  in  addition  to  conventional  boundary -layer  pressure 
measurements  for  experimentally  determining  friction  coefficients 
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in  flows  with  pressure  gradients .  The  disagreement  between  the 

experimental  and  empirical  friction  coefficients  wa3  probably  « 

caused  by  secondary  flows  in  the  boundary  layer. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  April  28,  1950. 
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APPMDIX  -  SYMBOLS 

The  following  symbols  are  used  in  this  report: 

local  coefficient  of  friction  based  on  density  evaluated 
at  stream  edge  of  boundary  layer,  T /p]_u-f 

,  local  coefficient  of  friction  based  on  density  evaluated 
at  -wall,  T/p  u-i 2 


(cf  ^x, 1 


(cf )©,1 


(°f 


empirical  local  coefficient  of  friction  based  on  Be£ 
and  density  evaluated  at  stream  edge  of  boundary  layer 


;f  )jc  w  empirical  local  coefficient  of  friction  based  on  Bej 


and  density  evaluated  at  wall  ’ 

empirical  local  coefficient  of  friction  based  on  BeQ  n 

t  y  j  Jl 

and  density  evaluated  at  stream  edge  of  boundary  layer 

empirical  local  coefficient  of  friction  based  on  Keg  v 
and  density  evaluated  at  wall  9 

turbulent-f orm.  parameter  assuming  variable  density  in 
boundary  layer ,  S*/e 

turbulent -f orm  parameter  assuming  constant  stream  density 
in  boundary  layer , 

Mach  number 

U  /vxl/W 

power -nr of lie  parameter.  —  =  / 'h\ 

u5  UJ 

total  pressure ,  absolute 
Prandtl  number 

static  pressure,  absolute 


Re5c.l 


Beynolds  number  based  on  equivalent  length  of  run  and  - 
kinematic  viscosity  evaluated  at  stream  edge  of  boundary 
layer 


Re?.w 


Reynolds  number  based  on  equivalent  length  of  run  and 
kinematic  viscosity  evaluated  at  wall 
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Ree,i 

Ee0/w 

r 

T 

u 

u5 

x 

X 

y 

7 
5 

8 

5* 


7lr 

0 


Reynolds  number  based  on  momentum  thickness  and  kinematic 
viscosity  evaluated  at  stream  edge  of  boundary  layer 

Reynolds  number  based  on  momentum  thickness  and  kine¬ 
matic  viscosity  evaluated  at  wall 

apparent  radius  of  radial  flow- 

temperature 

< 

velocity 
0.99  u^ 

distance  downstream  of  geometric  nozzle  throat 
equivalent  length  of  turbulent  run 
distance  normal  to  wall 
ratio  of  specific  heats 
boundary-layer  thickness 


boundary -layer  thickness  at  u§ 


displacement  thickness  assuming  variable  density  in 

05 


boundary  layer, 


displacement  thickness  assuming  constant  stream  density 


in  boundary  layer, 


temperature  recovery  factor 

momentum  thickness  assuming  variable  density  in  boundary 
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U 

P 

T 


momentum  thickness  assuming  constant  stream  density  in 


■boundary  layer, 


U 


kinematic  viscosity  of  air 
density 

shear  stress  at  vail 


Subscripts : 

0  value 
1  value 
v  value 


taken  in  surge  chamber 

taken  at  stream  edge  of  boundary  layer 

taken  at  vail 
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Figure  2*  -  Test— channel  instrumentation, 


Figure  3.  -  Pressure  probes. 
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(a)  Along  center  line  of  bottom  wall  of  channel. 
Figure  4.  -  Mach  number  distribution. 
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Figure  6*  -  Boundarylayer  profiles  on  side  and  bottom  walls  at  inlet  pressure  of  37  inches 

of  mercury* 
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Figure  6*  -  Concluded,  Boundary-layer  profiles  on  side  and  bottom  walls  at  inlet  pressure 

of  37  inches  of  mercury. 
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Figure  7*  -  Boundary-layer  development  along  channel 
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Distance  from  throat,  x,  in. 
(b)  Center  line  of  side  wall 
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Figure  8*  -  Concluded.  Turbulent -form-parameter  variation, 


channel  at  inlet  pressure  of  37  inches 


Turbulent  boundary- 
layer  development 

Theoretical 
(from  reference  3) 

- Using  (cf w 

-  -  Using  (cf )- 
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Figure  9.  -  Concluded.  Comparison  of  experimental  and  theoretical  boundary-layer  development  along  channel  at  inlet 

pressure  of  37  inches  of  mercury. 
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Figure  10.  -  Concluded#  Comparison  of  experimental  and  empirical  power-profile  parameters  along 

channel  at  inlet  pressure  of  37  Inches  of  mercury# 
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Figure  14.  -  Secondary  flows  in  side-wall  boundary  layers 
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